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I. INTRODUCTION 



Multi-component (or coupled) nonlinear PDEs had been a subject of considerable interest 
for many years (see, for instance, Ref. 1, and references therein). Recent revival of interest 
in multi-component PDEs is due to new discoveries and technological advances in nonlinear 
optics and physics of condensed matter. An important example is the incoherent spatial 
optical sohtons, or self-trapped spatially incoherent light beams, recently experimentally 
observed in nonlinear media [2], which are described by the multi-component nonlinear 
Schrodinger (NLS) equations [3]. Another example of possible apphcation of the coupled 
NLS equations is the creation and dynamics of solitary waves in the multispecies Bose- 
Einstein condensates [4] . Similar models of coupled nonlinear PDEs appear in the wavelength 
division multiplexing, i.e., copropagation of pulses in an optical fibre on beams with different 
wavelengthes [5-7] and in other important applications [8,9]. 

Some of the multi-component models are integrable. Integrable multi-component PDEs 
have another specific feature, which makes them important for applications as zero- 
approximation models for analytical description of the real phenomena. It has been known 
for quite some time that dimensional reductions of matrix generalizations of the integrable 
PDEs, such as the NLS and KdV equations, can produce a variety of new integrable equa- 
tions [10]. For instance, some of the coupled NLS equations are integrable reductions of 
the general matrix NLS equation. The A^-dimensional matrix NLS equation is the sim- 
plest integrable PDE associated with the {N + 1) -dimensional Zakharov-Shabat spectral 
problem [11-13]. Recently a variety of integrable coupled higher-order NLS equations was 
discovered [14-16], which are important in view of applications to the soliton propagation 
of sub-picosecond pulses in optical fibre [17-19]. Some of these integrable PDEs arise as 
dimensional reductions of the matrix complex modified KdV (cmKdV) equation, which is 
also associated with the Zakharov-Shabat spectral problem. 

In most cases, the multi-component PDE is not integrable. However, frequently the 
terms destroying integrability contain small parameters and the non-integrable equation can 
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be considered as a perturbation of the integrable one. In this case, a perturbation theory 
is required for analytical description of the effect of small perturbations. For instance, 
one is especially interested in the dynamics of solitons in nearly integrable PDEs. Soliton 
solutions to multi-component equations have many parameters, and their evolution may 
exhibit a variety of new interesting regimes. Therefore, it is necessary to have at hand a 
perturbation theory for multi-component nearly integrable equations. Such perturbation 
theory is developed in the present paper. 

Perturbation theory for nearly integrable PDEs has a long history [20-45] . There are two 
basic approaches in the perturbation theory based on the 1ST method. The first one origi- 
nated in works of Kaup [20] and Karpman and Maslov [21], where the perturbation theory 
was developed for nearly integrable PDEs associated in the integrable hmit with the 2x2 
matrix Zakharov-Shabat spectral problem. Quite different approach originated from works 
[22-26]. It was found that an integro-differential operator, generating the whole hierarchy 
of integrable PDEs related to a given spectral problem, called the recursion operator, has a 
complete set of eigenfunctions, which can be used for the perturbation expansion. Several 
other methods, not related to the 1ST, were apphed for description of the perturbed sohton 
dynamics. For instance, a method based on the Green functions was developed in Ref. 29. 
The IST-independent perturbation theories for solitons are usually referred to as the direct 
perturbation theories (see for instance, Ref. 40 and references therein). However, notwith- 
standing the long history of the perturbation theory, with rare exceptions, only the 2x2 
matrix spectral problems were considered. It was noted that construction of the perturba- 
tion theory for higher-dimensional matrix spectral problems along the standard approach 
becomes technically more involved. 

To overcome technical difficulties of the standard approach when dealing with multi- 
component PDEs, the method based on the Riemann-Hilbert (RH) problem was proposed 
in Ref. 41, where the perturbation theory was developed for the Zakharov-Shabat spectral 
problem of an arbitrary matrix dimension. The RH problem was used before for construction 
of the perturbation theory for the Landau- Lifschitz equation [30], the NLS and Maxwell- 
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Bloch equations [34], which are integrable by the 2x2 matrix spectral problem. The 
approach of Ref. 41 was applied to the Manakov system [42] (i.e., the two-component NLS 
equation), modified NLS equation [43] and massive Thirring model [44]. These examples 
demonstrate that the perturbation theory based on the RH problem always works. Recently, 
the RH problem was applied to nearly integrable equations on the half line, arising from 
the singular dispersion relations [45]. In Refs. 41-45 the perturbation- induced evolution 
equations for the spectral data were derived with the help of some matrix functional (below, 
the evolution functional Il{x,t,k)). It is important to emphasize that the form of the 
evolution functional is invariant under the gauge transformations of the considered PDE [43]. 
Thus, once constructed, the evolution functional is valid not only for the whole hierarchy 
of PDEs associated with a given spectral problem, but also for their images under the 
gauge transformations. Writing the dispersion law, generating the spectral problem, in an 
abstract form A(A;) (see the next section for details) we discover that the form of the evolution 
functional remains invariant under the change of the spectral problem as well. This invariance 
trivially extends to the general initial-boundary value problems. For instance, for the half- 
line, where one would expect a difference, we have found that the evolution functional has 
similar form [45]. Therefore, it seems that the approach based on the evolution functional 
is universal for construction of the perturbation theories for nearly integrable PDEs. It is 
also technically simple. Derivation of the perturbation-induced evolution equations for the 
spectral data using the evolution functional reduces to calculation of integrals. 

This paper is a further development of Ref. 41. The previous results are substantially 
advanced. In particular, the evolution equations for the spectral data are considerably 
simplified with the help of some identities found for the evolution functional. We start 
section y with a brief discussion of the multi-component integrable PDEs associated with 
the Zakharov-Shabat spectral problem. We consider two examples, the matrix NLS and 
cmKdV equations, however, our approach is valid for many other multi-component PDEs. 
We have not made an attempt to give a complete exposition of the properties of integrable 
equations. We need only the Lax representation. Hence a way of deriving the Lax pair for 
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an integrable PDE from the dispersion relation of its linearization is briefly indicated. For 
completeness of the exposition, a detailed derivation of the RH problem is given in section 
pn| . Solution of the RH problem for multi-component equations involves some technicalities, 
which are discussed and detailed derivations are provided in the Appendices. We derive 
evolution equations for the spectral data with account for perturbations in section |V]. There, 
for an important special case of the vector nearly integrable PDEs the main result of this 
paper is formulated in a theorem. In the final section the equations of the first-order 
perturbation theory for a single vector soliton are given in explicit form. 

II. PRELIMINARIES: INTEGRABLE MULTI-COMPONENT NONLINEAR 

PDES 

Here we briefly discuss integrable PDEs with emphasis on the multi-component equations 
whose reductions are important for applications. In particular, we consider the matrix 
nonlinear Schrodinger and complex modified Korteweg-de Vries equations. We do not try to 
review this subject, for general considerations the reader can consult, for instance, Refs. 1, 
10-13, 46-53 and the references therein. The purpose of this section is to remind some of 
the basic notions in the 1ST method. Though the approach below can be applied to any 
nonlinear PDE solvable by the RH problem, we restrict the consideration to the integrable 
equations associated with the A^-dimensional Zakharov-Shabat spectral problem ([l|). 

Consider the integrable PDEs which arise as the compatibility condition for the following 
N X N matrix linear system (Lax pair) 

9^.$ = ik[A, $] + iQ{x, t)$ = ^A{k) + U{x, t, k)^, (1) 
dt^ = iuj{k)^A + V{x, t, = m{k) + V{x, t, (2) 

with 
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where A(fc) = —ikA and Q = iuj{k)A are the dispersion laws, Q is called the potential. Here 
the overline does not mean complex conjugation by default, e.g., in general, the functions 
Qij and q^j are not considered as complex conjugate to each other. When the overline does 
denote complex conjugation in the text below, each time it will be specially indicated. This 
special case corresponds to the Hermitian potential Q, = Q or = q, and it will be 
referred to as the involution. The temporal evolution equation (^) is specified by choice 
of the dispersion relation uj{k) in the following manner. For simplicity, let the dispersion 
relation be polynomial uj{k) = Y^pLi uipk^, then 



M-l 



v{k) = -r{m^-^} = -n{k) + J2 Vpk"- 

Here the matrix function $(A;) is expanded into the asymptotic series: 



(3) 



k ^ 00, 



and the operator V takes the polynomial in k part of ^ on the asymptotics. For 

example, the Zakharov-Shabat spectral problem (P is derived in this way 



U = -P{$A$-^} = ikA + i[<l>'^^\ A] 



with the obvious relation 



g = [$«,A]. 



(4) 



Hence, the Lax pair satisfies the property TrU = — TrA and TrV = —TiQ. 

The integrable nonlinear PDE related to the Lax pair (|l])-(0) is given by the compatibility 
condition (in our case, polynomial in k) 



idtQ - d^V + [ikA + iQ, V] = 



(5) 



via setting k = 0, while the positive powers of k supply the expressions of the coefficients 
Vp in (^) as functions of the potential Q and its x-derivatives, Vp = Vp{Q,Qx,Qxx, ■■■)■ 
For instance, choosing uj{k) = we obtain the well-known matrix nonlinear Schrodinger 
equation. Indeed, in this case 

V = -2ik^A - 2ikQ - AQ^ + iAQ^ 



and the equation (B) becomes 



tAQt + Qxx + 2Q^ = 0. 



(6) 



For a complete classification of the matrix integrable NLS equations with various reductions 
to Hermitian symmetric spaces consult Ref. 10. A particular important case of equation 
(^) is the vector NLS equation, a generalization of the two-component vector NLS, which 
was shown to be integrable by Manakov [12]. The vector NLS equation corresponds to the 
Hermitian potential and the reduction n = N — 1 (see the expression for A). In this case 
we have 



^ . . . gi ^ 
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and matrix equation (P) becomes the vector NLS equation: 
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(7) 



tdtqi + dlqi + 2{Y,\q^\^]qi = 0, 1 = 1, 



,n. 



(8) 



Let us consider another important example of mult i- component nonlinear integrable 
equations. It is given by setting uj{k) = 4k^. After simple computation we get 

V = -Ate A + Aik^Q + 2ikA (iQx + Q') - iQ^x + [Q, QA - 2iQ\ 



which produces the matrix cmKdV equation 



Qt + Qx.. + 3 (q^Qx + g.Q') = 0. 
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(9) 



Combining together the considered two dispersion laws, i.e., letting uj{k) = 4iek^ + 2if3k'^, 
one can obtain 

iQt + (3 A (q^^ + 2g3) + ieQ:,^^ + 2>ie {Q^Q^ + Q^Q"^) = 0, 

a special case of the (generally, non-integrable) matrix higher-order NLS equation 

iE, + A{aiErr + a2^^) + t [a^Errr + a4(^^)r + a^{E'^)rE] = 0, (10) 

Our examples are just illustrative. There is many other integrable matrix PDEs which 
we do not mention here. However, the perturbation theory developed in section |^ applies 
to such PDEs also. 

III. RIEMANN-HILBERT PROBLEM FOR MULTI-COMPONENT PDES 

In this section we derive the RH problem for the multi-component integrable PDEs and 
discuss its solution (for more details consult Refs. 46, 54-57). We are interested in the initial- 
value problem for nonlinear PDEs on the whole real line with the asymptotically vanishing 
conditions (Cauchy problem) qij ^ as |x| oo. The vanishing asymptotics allows us to 
concentrate entirely on the spectral equation (|l]) in derivation of the RH problem, while t- 
dependence enters parametrically in our approach (for the RH problem for initial-boundary 
value problems consult Ref. 47). Below we omit the explicit t-dependence for simplicity of 
the presentation. To begin with, let us summarize the properties of the spectral problem 
(0). Define the following N x N matrix projectors 

Hi = diag(4, 0), H2 = diag(0, In-u). (11) 

Then A = Hi — H2 and any matrix can be decomposed into the sum of two matrices, 
commuting and anti-commuting with A: 

$ = $W + $W, ^^""^ = Hi^Hi + H2^H2, ^^"^ = Hi^H2 + H2<^Hi, (12) 
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where = and {v4,$('^)} = 0. We will use the block-index notations for the 

decomposition of matrix $ with respect to the projectors Hi and H2: 



The RH problem is a one-to-one mapping (nonlinear Fourier transform) between the 
set of smooth (e.g., belonging to the Schwartz space) potentials Q{x) and some set of the 
spectral data. To identify the RH problem one must construct two solutions, one, $+(a;, k), 
to the spectral equation ([l|) and the other, $l^(x, fc), to the adjoint equation (the second 
function is inverse of some matrix function satisfying ([|), hence the in its definition), 
holomorphic with respect to the spectral parameter k in some complementary domains 
covering the whole complex /c-plane. Such solutions can be built from the columns and rows 
of the Jost solutions J±, i.e., solutions defined by the asymptotic conditions: J±{x, k) —>■ I 
as X ±00. 

As TtQ = 0, letting x ±00 we conclude that det J-t = 1. Hence, the columns of either 
of the two Jost solutions give a linear basis in the space of solutions of the spectral equation. 
The inverse matrices J^^ satisfy the adjoint spectral equation: 

d^^ = ik[A,^]-i^Q. (13) 

Each column of J±{x, k) and, respectively, row of J^^{x, k) is holomorphic and bounded in 
either upper (ImA; > 0) or lower (ImA; < 0) half of the complex /c-plane. Indeed, this can 
be easily seen from the Volterra integral equations for the Jost matrices written for the two 
blocks of columns: 

X 

Mx,k)Hi = Hi+t J dee-2^'=(^-«)^^Q(0J±(e,A;)i/i, 
±00 

X 

Mx,k)H2 = H2 + t J dee2^^(^-^)^^Q(0J±(e,A;)i/2. 

±00 
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The columns of J+{k)Hi and J-{k)H2 are holomorphic and bounded in the upper half of 
the complex plane, while columns of J^{k)H2 and J_{k)Hi have the same property in the 
lower half plane. Similarly, the rows of J^^ satisfy integral equations 



X 

H^JZ\x,k) = H^-i J d^H^Jz\^,k)Q{Oe 

±00 



X 

H2Jz\x,k) = H2-i J d^H2Jg\^,k)Q{0e-'''^'-^^''\ 

±00 

from which we immediately conclude that HiJ^^{k) and H2Jz^{k) are holomorphic and 
bounded in the upper half plane, while H2J^^{k) and HiJz^{k) have the same properties 
in the lower half plane. 

On the real line, the Jost solutions are transformed into each other by the scattering 
matrix S{k), 

J_(x, k) = J+(x, A;)e^'=^^5(A;)e-^'=^^. (14) 

For the Hermitian potential, Q'^ = Q, the matrix Jost solutions satisfy the involution 
(here the overline means complex conjugation) 

Ji{x,k) = J^\x,k), (15) 

where the spectral parameter takes complex values in the upper or lower half plane depending 
on the considered column of the Jost matrix. In this case, the scattering matrix also satisfies 
the involution 

^t(^) = S'\k), k G Re. (16) 

The holomorphic matrix functions ^+{k) and $I^(A;), satisfying equation (0) and (|13|), 
respectively, are given in terms of columns and rows of the Jost solutions: 

$+ = J+Hi + J^H2, = HiJ-^ + H2Jl\ (17) 
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The above defined matrix functions are liolomorpliic and bounded in tlie upper and lower 
lialf planes, respectively. They have the following asymptotics 

^±{x,k)^I, k^oo, (18) 

which follow from the Volterra integral equations for the Jost solutions. For the involution 



(|15|) the matrices ^+{k) and ^_^{k) are related via 

<l>l{k) = <^z\k). 



(19) 



These matrices can be conveniently expressed in terms of only one Jost solution and elements 
of the scattering matrix S{k). Indeed, 



$+ = J+e"'''^{Hi + SH2)e 



-ikxA 



ikxA 



Denote Sj^ = Hi + SH2, = H2 + S^^Hi. These matrices provide a factorization of the 
scattering matrix: = SS^. Similarly, S+ = Hi + H2S~^ and S- = H2 + HiS, which also 
factorize the scattering matrix: S^S = Then 



ikxA o „—ikxA 



ikxA o „—ikxA 



-1 ^ikxA^ ^—ikxAj—1 _ ^ikxA ^—ikxAj—l 



(20a) 
(20b) 



The factorization matrices have the block-triangular structure. For instance, 5"+ and 
are upper and lower block-triangular, respectively: 



















^0 aj 









(21) 



where b = 5'/^/, a = Suji, h = {S ^)nj, and a = (S* ^)n,ii- The following identity follows 
from these definitions 
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bb + aa = In-u, k G Re. (22) 

For the involution (0) the factorizations satisfy 

Si{k) = S±{k), ke Re. 

Hence, b = b^^ and a = aMn the case of involution. 

Considering the product we obtain the problem of analytic factorization of a 

matrix G{k) given on the real line, i.e., the matrix RH problem: 

A;)$+(x, k) = e''="^G'(A;)e-'^"^, A; G Re (23) 

and A;) / for A; oo. Here the matrix G = S+S+ = S^S^ reads 



G 



<i. b ^ 

y b /jv-n J 



(24) 



As it was mentioned, the RH problem is a nonlinear mapping between the potential 
Q{x) and the set of the spectral data, which are necessary for unique identification of the 
solution to (|23|) . For instance, given a potential, one can obtain the spectral data by solving 
the spectral equation and its adjoint for $+(x, k) and $I^(x, A;). Conversely, by asymptotic 
expansion of ^±{x,k) as A; ^ oo one recovers the potential. The asymptotic expansion of 
^^(x, k) can be derived via integration by parts (in the blocks with e='=^*'^^) in the Volterra 
integral equations for J± and J^^. It reads 

(X oo \ 

— oo X / 

Hence, we obtain (cf. with (^) 

Q{x) = lim k[^+{x,k),A] = lim k[A,^z\x,k)]. (25) 

k^oo k—*oo 
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Solution of the RH problem: Normalization 



Let us discuss the way of solution of the RH problem. The coordinate dependence is not 
important for this purpose and it is omitted below. In general, the determinants det ^+{k) 
and det have zeros and the RH problem is said to be non- regular or with zeros. Note 

that the determinants do not depend on x as it readily seen from equations (|T]) and (P!3[). 
We consider only the RH problem with zero index, i.e., 

oo 

J dln{detG'(A;)} = 0, 

— CO 

assuming that G{k) is non-degenerate. Since <l>± / as — > oo, the index is equal 
to the difference between the number of zeros in the upper and lower half planes. Let 
ki,...,kM and ki,...,kM be zeros (where some may be equal) of det<l'+(A;) and det^Z^{k), 
respectively. Two conditions are imposed on the zeros. First, the geometric multiplicity of 
a zero must be equal to its order (which we will refer to as the algebraic multiplicity). Here 
the geometric multiplicity of kj is defined as the dimension dj of the null space of $+(/cj), 
i.e., dj = N — rank$_|_(A;j). Trivially, the two multiplicities are equal to 1 in the case of 
simple zeros. In the case of involution zeros of the RH problem come in complex conjugate 
pairs (due to formula (|19])), with coinciding algebraic (and geometric) multiplicities within 
each pair. (Due to the involution the index is equal to zero, in this case detG'(A;) is real 
and definite.) Second, as we are mainly interested in the Hermitian potentials Q, i.e., in the 
case of the involution, we will consider only paired zeros kj and kj, j = 1, s < M, whose 
algebraic multiplicities are equal, however without assuming them to be complex conjugate 
to each other. The algebraic multiplicity Uj of kj always satisfies the following inequality 
(see Appendix A) 

Vj > N — rank$+(A;j) = dj, 

while rank$+ > max(n, N — n) by construction of $+ (|T7|). For instance, if n = A^ — 1, there 
can be not more than one vector in the null space of ^^(A;), i.e., dj = 1 for all j = 1, s. 
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Hence, for this reduction, zeros of det ^+{k) must be simple to satisfy the equal multiplicity 
condition. Similar for zeros of (More detailed consideration of the multiplicities of 

zeros is placed in Appendix A.) 

Let det^+{k) and det$I^(/c) have zeros ki,...,ks and ki,...ks, respectively, with the 
(algebraic and geometric) multiplicities Exactly Vj independent vector-columns 

\Pi^) (vector-rows (pp^l), where / = 1,...,^'^, correspond to each zero kj (respectively, /cj), 
such that 

*+(%)IpP) = o, (#|$:^(^,) = o. (26) 

To specify a unique solution to the RH problem, additionally to the continuous datum G(/c), 
the set of the discrete data {/c^, (pp^M — ■■■■i^j-ij — ^i -i^} must be given. This 

becomes evident from the fact that the RH problem with zeros (i.e., non- regular) can be 
regularized, i.e., reduced to a modified RH problem without zeros by factoring them out 
with some rational matrix T{k): 

$+(fc) = 0+(fc)r(fc), ^Z\k) = T-\k)r-}{k). (27) 

The regularized RH problem reads 

(j)Z^{k)(j)+{k) = T{k)e'^''^G{k)e-'^''^T-^{k), A; e Re (28) 

and 4>±{k) ^ / as /c — > oo. 

The rational matrix functions r(A;) and r~^{k), (for details see Appendix B) 

^ A~T^ ' -^+2- -TTT. ' ^^^^ 

j,l;i,m '^J j,l;i,m '^J 

where 

rvj Ki 

have the same zeros as ^+{k) and $I^(/c), respectively, and the same null spaces: 

r(A;,)|pp^) = o, (#|r-i(^,) = 0, / = j = 

14 



Moreover, det T = 11^=1 (^^zf^^ • Precisely these properties allow us to factor out zeros of 
^+{x, k) and $I^(a;, k). 

It is important to emphasize that there is a freedom in choice of the basis vectors spaning 
the null spaces. Indeed, it is easy to verify that the regularization matrix F is invariant under 
the transformations: 

IpF^) - E \l^)Mi^. {v?\ - E M^lm. (30) 

m=l m=l 

where M^^^ and M^^^ are arbitrary non-degenerate uj x z/j-matrices. For the involution, due 
to equation (|T^) we can choose = bPV- This evidently leads to 

Ft(A;) = T-\k). (31) 

The (x, t)-dependence of the RH data can be found in the following way. The t- 
dependence of the continuous datum (which does not depend on x) can be derived from 
equations (p3D and (||). We obtain 



dtG=[G,n]. (32) 

This equation can be cast in a more explicit form: 

dth = -2iuj{k)h, dth = 2iuj{k)h. (33) 

Recalling that TiQ = and TiV = — Trfi, one can verify that the determinants of $± do 
not depend on the co-ordinates. Therefore, the zeros are co-ordinate independent as well: 

d^kj = dtkj = 0, d^^j = dtkj = 0, j = 1, s. (34) 

Now, let us derive the co-ordinate dependence of the vector-parameters. Differentiating 
equation (|26| ) we obtain 

- ^k,A\p'f'>)) = 0, {d.lp'f^) + n{kj)\p'f^)) = 0. 

Therefore 

15 



m=l m=l 

where F and K are some matrices. Using a suitable invariance transformation of the type 
(pop one can put F = K = Q without loss of generality. Hence we have 

d^\p?) = ^k,A\p'f^), d,\p'f^) = -fi(A:,)bFV (35) 

Similarly 

dAP?^\ = -{P?^\^k,A, = {p^^\n(k,). (36) 

Some comments are needed on the reconstruction of the potential Q. First of all, the 
soliton part of the potential is given by the rational matrix function T{k) (p9|). The pure 
soliton potentials are also called reflectionless, they solve the simplest RH problem with 
G = I, i.e., with zero reflection coefficients: b(A;) = b(A;) = 0. The discrete RH data 
have the following meaning. Zeros provide the amplitudes and velocities of the solitons, 
while the null vectors give their initial position, polarization and phase parameters. Second, 
the radiation part of the potential is given by the solution to the regularized RH problem 
The continuous RH data h{k) and h{k) represent the nonlinear spectral densities 
of radiation. The regular RH problem is equivalent to some matrix integral equation of 
the Fredholm type (for instance, consult Ref. 46). Though the solution cannot be given 
in explicit form, its properties can be explored by the standard technique of the theory of 
Fredholm integral equations. 

For the following, it is convenient to introduce x-independent null vectors l-P/"'^) and 
(P?^| setting 

Then, we have an x-independent set of the RH data, which we will call the spectral data: 
{h{k)Mk),kj,\Pl'),kj,(Pi'''\,l = l,...,u,,j = l,...,s}. 

As the illustrative example, consider one-soliton solution. It is given by r(A;) having only 
one pole, say ki = irj + C,- In this case 
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Lm=l 



with Aim = {Pi\Pm)- Here {pi\^_^{ki) = and $+(A;i)|j9;) = 0. In the particular case of 
n = N — 1 and the Hermitian potential we get the vector soliton solution, for which 

Pr = f^, \p)=exp{zhxA}\P). 

Define complex parameters Ci = Pi/ Pn- The t-dependence of Ci then follows from equation 



(^). It can be accounted for in a convenient way by introducing real t-dependent parameters 
X, and 5i, I = 1, n, the soliton position and phases of its components. Let Ci = Oie^^^'^^^^ , 
where 9i = sie^^K The amplitudes s;, satisfying 

n 

Y.sl = l. (38) 

1=1 

describe polarization of the multi-component soliton. Integrating equation ( P3| ) we obtain 
X = xo + 5l = 5l, + -\m{{^-^r^)u:{^r^ + 0}t. (39) 

7] 7] 

The vector soliton then takes the form 

qi = 2ir/^ze^(«/^)"sechz (40) 

with z = 2rj{x—x). For instance, for the vector NLS equation uj{k) = 2k'^ and we arrive at the 
n-component generalization of the well-known vector soliton solution for the two-component 
NLS equation [12]. 

IV. PERTURBATION-INDUCED EVOLUTION OF THE SPECTRAL DATA 

A perturbation of the integrable PDE, following from the Lax representation can 
be written in the form 

iQt -V, + [ikA + iQ, V] = R, (41) 
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where R (which is fc-independent in case of the Zakharov-Shabat spectral problem) represents 
the terms destroying integrablity. The perturbation matrix R, similar to the potential Q, 
satisfies [A, R] = 0. For the Hermitian potential Q'^ = Q we get also R^ = —R. 

Below we derive evolution equations for the spectral data with account for arbitrary 
perturbation. For simplicity of the presentation we will frequently omit the explicit depen- 
dence on X and t. To distinguish between the "integrable" and "perturbation" contributions 
to the evolution let us assign the variational derivatives to the latter. For instance, the 



perturbation-induced evolution in equation (^iD reads 



^5Q _ explicitly i-^ = r and i-r^ = — r, 

ot ot ot 

where r = Rjjj and r = —Rn,i- For Hermitian Q, r = 

We start with derivation of evolution equations for $_|. and $1^, the solution to the RH 
problem (^). In other terms, we are going to derive the generalized Lax pair (see equations 
(^) and (^) below) for the perturbed nonlinear PDE (0). Differentiation of equation (|l]) 
with respect to t gives 

="=i'^-^ri+^'««+^*- 



and consequently 



Let us integrate the above formula for $ = J±. We get 

Here we have used that 6J± — as x — > ±oo. Let us employ the relation between J± and 



to rewrite this formula in a more convenient form. Using ( |20a|) we get 

A T 

(x, k) = J±(x, A;)e*'=^^5±(A;)T(±oo, x; k)S^\k)e-'''''^. (42) 



6t 

In much the same way, using (|20b|) , we derive 



A 7-1 _ _ _ 

± ^x,k) = -e'^'''^Sl\k)T{±oc,x;k)S±{k)e~"'''^Jg\x,k). (43) 



5t 

Here we have introduced the notations: 

X 

T{±oo,x;k)= j Aie-'^^^^-^m+e'^^^, (44) 

X 

T(±oo,x;A;)= j e-'^^'^^Z^ m .e'''^'^ . (45) 

±oo 

These matrix functionals will enter every formula describing perturbation-induced evolution 
of the spectral data. For the involution, due to formula (0) these matrix functionals satisfy 

Tt(±oo, x; k) = -T(±oo, x, k). (46) 

Using equations (^2p, (^3]), and the definition of the scattering matrix S ( JT^ j we get 

_ _ e-''^"^- (J+ V_) e''^"^ = -5;iT(oo,a;)5+5 + SS-T{-oo,x)SZ^ 

= 5'+T(— oo, x)S'r"'" + S^^T(x, oo)S^. 

Setting X ±oo in this formula produces two simple equivalent formulae: 

= 5+(A;)T(-oo, oo; k)Sz\k) = S-\k)T{-oo, oo; k)S4k). (47) 

Now we can easily obtain the perturbation-induced evolution of the solution to the RH 
problem. Taking into account the definitions of S± and S± and (|20|) we write: 



^^gifcxA (_T(x, 00) + T(-oo, oo)i72) e-'^""^ = ^^g^'^^^ne"^'^^^, (48a) 



AkxA'Q -ikxA^'^+ I AkxA ^ / tt i rr c-l^ T-l^-ifczA 



gifcxA i^y^^^ _ /72T(-oo, 00)) $:ie-^^^'^ = -e'^^Aj^^-i^-ikxA_ ^^g^^ 
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The r.h.s.'s of the above formulae contain the evolution functionals, 

n(x, k) = — T(x, oo; k)Hi + T{—oo, x; k)H2, 



(49) 



n(a;, k) = HiT{x, oo; k) - /J2T(-oo, x; k), (50) 

which account for the perturbation-induced evolution of the solution to the RH problem. 
As follows from formula (^61) , for the case of the involution the evolution functionals satisfy 

U^k) = -U(k). (51) 

From the definitions (^), (|50|) and also (^If), (^) it is easy to see that the matrices 
^tkxAYi^j^^^-tkxA g^j^^ ^tkxAYi(^j^^^-ikxA ^^^^ meromorphic and bounded in the upper and lower 

half planes of the /c-plane, respectively. They have simple poles respectively at zeros of 
det$+(A;) and det$I^(A;). 

Let us write down the generalized Lax representation for the perturbed PDE (^T]) in 



terms of ^+{x, k): 

d,^+=ik[A,^+]+iQ^+, (52) 

dt^+ = + V<^+ + <^^e'''^Ajj^-i''xA_ ^53^ 

It is easy to check by direct calculation that the compatibility condition for the above linear 
system is equivalent to the perturbed equation (^) (indeed, d^H = e~'^''^^^^^ R^+e^'^^^) . 



Now, let us derive the perturbation-induced evolution of the spectral data. From (^Hj) 
we immediately obtain 

Hence, the complete evolution of the continuous datum reads 

dtG = [G, Q]+GU- UG, (54) 
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where we have taken into account the integrable evolution given by (|32D. It is important 
to notice that the l.h.s. of (|5^ ) does not depend on x. Therefore, we can put x ±00 in 
this equation to simphfy it. Equation ( |5^ ) can be rewritten in a more exphcit form. Using 
the definition (p^, formulae (^) and we get two equivalent evolution equations for b 
corresponding to the two limits x — ^ ±00: 

dth = -2iujik)h + bif2T(-cx), 00; k)H2 + i/iT(-cx), 00; k)H2, (55a) 

^tb = -2iuj{k)h - HiT{-oo, 00; k)Hih - ifiT(-oo, 00; k)H2. (55b) 

Similarly, 

dth = 2iu{k)h + bi7iT(-oo, 00; k)Hi + i/2T(-oo, 00; k)Hi, (56a) 



dth = 2iuj{k)h - //2T(-oo, 00; k)H2h - if2T(-oo, 00; k)Hi. (56b) 

Evolution of zeros kj and kj of det$+(/c) and det$l^(/c), respectively, is derived by 
differentiation of the determinants. In the integrable limit zeros do not depend on t. Then, 
for instance. 



dkj _ dtdet<^+{k) 



dt dkdet^+{k) 



Tr{n(A;)}det$+(A;) 
~ dkdet <l>+{k) 



To calculate the r.h.s. in this formula recall that the evolution functional 11 has simple pole 

ai k = kj and that det $+(A;) can be written as 

det <I>+(A;) = det 0+(A;) det T{k) = det 0+(A;) [] 7 — ^ , 

i=i\k — kij 

where det0+(fc) 7^ 0. Simple calculations give 

dkj _ TT{ResU{kj)} 
dt Vj 

Here "Res" denotes the residue of n(A;) at A; = kj. Noticing that the l.h.s. does not depend 
on X, we simplify the above equation setting x — » ±00: 
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dkj Tr{ResT(— oo, oo; /Cj)if2} _ Tt {ResT {— 00,00; kj) Hi} , . 

dt Uj Uj 

Similarly, using det$I^, we get 

dkj Tr{ResT(— 00, 00; /Cj)if2} _ Tr{ResT(— 00, 00; A;j)i7i} 



dt Uj Uj 



(58) 



In derivation of the perturbation-induced evolution of the null vectors we will use the 
following remarkable identities (written for the x-independent null vectors, see (P?])) 

Resn(fc,)|P/^'^) = -^1^/^'^), (PF^|Resn(fc,) = (Pp)|^, (59) 

as well as other two identities: 

$+(%)e*^^^'^Resn(A;j) = 0, Resn(kj)e-''^^^Hz\kj) = 0, (60) 

which follow from (^81). To verify the identities (^) one can proceed as follows. Introduce 
functions Fj^\k) = {k ~ kj)^:^^{k) and F-'\k) = {k - kj)^_{k). The regularization 
matrices r~^(A;) and T{k) have simple poles at k = kj and k = kj, respectively. Hence the 
introduced matrix functions are holomorphic in some neighborhoods of these points. Now 
compute the product 



= e"'^''^ResU{kj)e-"'^''^\pY>). 

On the other hand. 



St St +'■''■ 

Multiplication by of the latter formula and comparison with the former leads to the 
first identity in (^). The second one can be checked in similar way. 
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Now let us derive evolution equations for the x-independent null vectors defined in 
To this goal we simply differentiate equation (p6D 



Denote ^r{kj) the value of the regular part of n(fc) at = kj, 



U^ik,) = \ Uik) - ^ . (61) 

'^j ) k=k. 



Then, using (|59D and (^) to cancel out secular terms, we arrive at 

^^{k,)e^'^-4^-^ + n.(A:,)|P/^M = 0. 



Using the same arguments as in section |ITl|for derivation of the integrable evolution we get, 
adding the latter, 

d|P/^-)) 



dt 



Similarly, 



'^^^'^'^'^ -{p?m,) + {p^p\ii^{k,), (63) 



di 



where 



Resn(A;j^ 



n,(A:,) = jn(A:)- _ • (64) 

In the case of the involution equation (|63|) is Hermitian conjugate to (|62|) . Note that the 
l.h.s.'s of equations (^) and (|6^) do not depend on x. Hence we can sent x ±oo to 
considerably simplify these equations. Consider, for instance, (|62D. Letting x ^ oo and 
introducing the notations 
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and T^^) = H1TH2, T^"^^ = H2TH2 we get the following system 



^-L^ = -^u;{k,)\P,^^\ 1) - T«(-oo, 00; k,)\Pp\ 2), 



{za;(A;,-)-T(2)(-oo,oo;fc,)}|P/^\2). 



(65a) 



(65b) 



Here Tr{kj) is value of the regular part of the matrix T{k) at k = kj-. 

^ ResT(A;j) 1 



T,(A;,)= T(A:)- 



(66) 



When dealing with vector PDEs, i.e., for n = N — 1, one can define the polarization-phase 
parameters of vector solitons as quotients of components of the null vectors (note that the 
zeros are simple in this particular case). Let C^''"'^ = P^''^ /P^\ I = 1, ...,n = N — 1, where 
|p(i)) = (pI^\ ...,pf^^f. Then from (H) we obtain: 

dC/^) 



dt 



{—2iuj{kj) + Tr7VAf(— 00, 00; kj)} Ci^^ — TriN{—oo, 00; kj), I = 1, n. 



This particular case {n = N — 1) contains the vector NLS (|]) and the complex modified 
KdV (1^) equations as examples. In view of considerable importance of such vector nonlinear 
PDEs, we formulate the result of this section in the following theorem. 



Theorem: Let 



iQt ~l~ ^oi^Q y Qxi Qxxi •••) 6-R(x, t, Q, Qxi Qxxi •••)' 



(67) 



he a perturbed nonlinear PDE associated with the N x N matrix Zakharov-Shabat spectral 
problem, 



d^^ = ik[A, $] + iQ{x, t)^, A = diag(4, -1) 



(68) 



where (n = N — 1) 

Q 



(0 




( r) 








, q 
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Here the matrix function Vq represents the limiting integrable evolution given by the dis- 
persion relation Q{k) = iuj{k)A, while R contains the terms destroying integrability (e is a 
small parameter) . Then the perturbed evolution is equivalent to the following evolution of 
the spectral data: 



Here 



dt 



—eResTiy]\r{t, kj), 



— ^ = {-2iuj{kj) + eTrNNit, kj)} Cj;^'' - eTriN{t, kj), I = 1, n, 



dbijk) 
dt 



{-2iLj{k) + eT^N^t, k)} bi{k) + eT,Ar(t, k), / = 1, n. 



T{t,k)= / dxe-"'^'^^~'{x,t,k)R{x,t)^+{x,t,k)e 



ikxA 



(69) 



(70) 



(71) 



(72) 



T,(t,A;,) = \Tit,k) 



ResT(t, kj) 



k-k 



The matrix function $_|_(x, t, /c) solves the spectral problem and the RH probl 



lem 



^\{x, t, k)^+{x, t, k) = e*^^^G'(A;, t)e-*'^^^, k G Re, 



-ikxA 



^+{k) ^ I, k ^ oo 



of analytic factorization of matrix G{k, t), 

[ 

G{k,t) = 



V 



b = (61,62, ■■■,bn) 



I 



In b(A;,t) 
h\k,t) 1 

The Riemann-Hilbert problem has (simple) zeros kj{t), j = l,...,s, given by a{kj,t) = 
det^+{x,t,kj) = 0. The vector- columns \C^^\t)) = (c[^\t),dj\t), ...,Ci,^\t),l) are the 
null vectors of the matrices $+(x, t, fcj)e 



$+(a;,t,A;,)e^'=^^^|C(^)(i)) =0, J = l,...,s. 



25 



The initial spectral data are obtained via solution of the spectral equation ( [g^ j and represent 
the spectral characterization of the potential Q{x,0). For real k, the spectral densities of 
radiation bi{k,t) and the function a{k,t) satisfy the following identity 

n 

aa = 1 — y^^bibi. 

1=1 

The potential q{x, t) is reconstructed by the formula 

qi{x, t) = 2 lim k{^+)i^N{x, t,k), I = l,...,n. 

fe^oo 

Some comments are necessary on the use of equations (|69D-([n]) for the spectral data. 
These equations are exact, i.e., they account for the perturbation exactly. As a consequence, 
for the non-integrable PDE (0), these equations are non-closed ODEs. Equations (|69|)-(|7lD 
are non-closed because they contain explicitly the matrix $+(a;, t, fc), solution of the RH 
problem, obtaining which requires knowledge of the spectral data. Therefore, equations 
(|69|)-([7T|) serve as the generating equations for the perturbation expansion: expanding the 
spectral data into the asymptotic power series in e, one obtains the sequence of closed 
approximate ODEs for the spectral data. On this way, one does not need to solve the 
RH problem - the computations are algebraic. In the next section we consider a single 
multi-component soliton as an example. 

V. MULTI-COMPONENT SOLITON UNDER PERTURBATIONS 

In this section we apply the theorem for construction of the perturbation theory for a 
single multi-component soliton. This can be done without specifying the dispersion rela- 
tion uj{k) determining the evolution of the spectral data in the unperturbed PDE. Hence 
our results apply to all nearly integrable vector PDEs associated with the Zakharov-Shabat 
spectral problem. Below we derive equations describing the evolution of the soliton param- 
eters and give formulae for the first-order radiation. For a single vector soliton given by 
formula (^Op, i.e., qi = 2ir]9ie'^^^^'^^^ sechz , where z = 2r](x — x), the regularization matrix F 
has the form 
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Here ki = ^+irj and we have used the vector notations \6) = (6'i, 9n)^ and {9\ = {9i, ...On)- 
To simphfy some of the calculations introduce the following basis (in the n- dimensional 
subspace) 



The explicit form of the vectors l^'''-') for / = 2, ...,n will not be needed at all. Also we will 
use the basis 



|ei) = (1,0,..., Of, lea) = (0, 1, 0, 0)^ 



\en) = (0,...,0,1)^ 



With the help of the unitary transformation matrix S, defined as 

n 



(74) 



1=1 



the regularization matrix F can be considerably simplified: 

/ 

f = ETE-^ = I 



17] 



{k — ki) cosh 2; 



e-"|ei)(ei| e^^^/^^^ld) 



(75) 



Note the evident property = T'^ . 

The transformation matrix E depends on t due to the perturbation-induced time depen- 
dence of the 6'-parameters, but does not depend on x. This simple fact allows us to use the 
transformation with E inside the integrals defining T (^). Taking into account that = F 
for the pure soliton solution of the unperturbed PDE, in the first order we get 



Here we have defined 



R = ERE 



dxe-*'=^^F-^i?Fe*'=^^. 



^ B\r)^ 



-(r|5t 
27 



(76) 



where, for convenience of the presentation below, we have changed the notation for the 
perturbation: |r) = r = (ri, ...,r„)^. 

In fact, we will need only one diagonal element 



' NN 



t™=(e-th)^^= / dx(ft/if)_ 

— oo 

and the following non-diagonal matrix elements 

n n 

'^IN = ("~^^")ijv ^ ^ ^fm^rnN = Ol^lN + BfJ^f^N = OifiN + Fi, 

m=l m=2 

where / = 1, ...,n. The second term simplifies as follows 



m=2 



m=2 



dxe-^^^^r^r, {ri-9i{9\r)). 



After simple calculations we get: 



^NN = T / d2;sech^^ 
4 J 



k — ki k — ki 



{ro{z)+ro{-z)) , 



(77) 



(78) 



+ 



e-2^(A; - k^) + e^^(A; - ' r^{z) \ + F^. 



J2z 



Here we have used the notation 



1=1 



(79) 



(80) 
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A. Evolution of the soliton parameters 



Let us first derive evolution equations for the soliton parameters rj, ^, x, and 6i, I = 
1, n. Using the identity —2ikiX+i{^/r])z = z+2r]x—2i^x and the definition Ci = ^^^e^*-''"*^^^ 
from section |T|, we obtain: 



oo 



ResTjv7v(^i) = ^ J dze ^sech^z (ro(-2) + ro(— 2:)) , (81) 



-00 



where 



1 °° 

—00 



TriMih) = Q {-2ixRes{Tr,N{ki)} + T,^^(fci) + Jq} + fi, (83) 



00 
1 r 

Jo = — J dzsech^^; (^coshz + ze~^^ (^0(2:) — ^o("~^)) 



and 



/; = F/(A;i) = 1 dzsech^ (e'^^/^N - ^^ro) 



4?7 

Let us first consider the more involved derivation of equations for x and 61. From equa- 
tions dZO), dH), and (H) we get: 

dCi 



dt 

from which it follows that 
d\Q\' 



-2iu{ki) + 2eixResTNN{ki) - eJo) Q - efi, (84) 



4Im{u;(A;i)} - 4exIm{ResT^^(A;i)} - 2eRe{Jo} \Ci\' - 2eRe{//C;} 



dt 

Recalling that ^27=1 = 1 we obtain: 



dt ~ 4r] ^~dt r/ dt' 
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dt V dt dt ) 

The rest calculations are straightforward substitutions and using the identity Ya=i fi^i — 0; 
which follows from the definitions of Q, tq, and After simple calculations one gets a 
system of equations for the soliton parameters: 



df] e 
dt ~ ~2 



oo 



dz sechzReIro}, (85) 



-oo 



dt 2 



oo 



/d.Bech.-tanhdm{r„}, (86) 



dOi ^ [ 2Im{A;ia;(A;i)} e /" , , -r^ r / , \t r t 

-— = zOi < / d2;sechz 42;Re|ro| + ri[l — 2;tanh2;)Im|ro| 

dt r] zr]^ J L 

^. — oo 

oo 

+ — J dzsechz {Oiro - e-^^^/^^V,) . (88) 

^ —oo 

These equations can be compared with the adiabatic equations derived by Karpman [21] 
for the single scalar soliton. First, it is convenient to introduce the average phase 6 of the 
soliton by setting 

1=1 

The evolution equation for the average phase then follows from equation (p8|): 

d5 2Im{A;ic<j(A;i)} e /" , i -r. r i i \t r t1 /or^^ 

— / dz sech^; 42;Ke|ro| + ?7(1 — 2;tanhz)lm|ro} . (89) 



oo 



-,2 

— oo 



dt 7] 27/2 

Remarkably, the slow evolution of the soliton amplitude r], phase gradient ^, position x and 
average phase 5 is given by equations similar to those derived for the single scalar soliton. 
The only trace of the vector nature of the soliton in equations (P5|)-(|S7D and (^Uj) is that the 
"scalar" perturbation tq obtains by averaging the original vector perturbation as follows 
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1=1 

The equation for 9i = sie''^' can be cast in the form of two separate equations, one for 
the polarization parameters si and the other for the phases Si. We get: 



dsi _ e 
dt At] 



J dzsechzRe{siro-e~'''^^'^>-'^'ri} , (90) 



d6i _ 2lm{kiuj{ki)} e 
dt T] 2rf 



J dzsech^; ^zRe{ro} + ^^(l — 2;tanh2;)Im{ro} 



-OO 



+ — J dzsech2lm{ro-sr^e-*(«/'')^-^'''n}. (91) 

^ — OO 

Note that the equation for 6i contains si in the denominator as a reflection of the fact that 
the phase is not defined for the components which are not excited. 



B. First-order radiation 



Now let us consider the evolution of the spectral densities bi{k) of radiation. Taking into 
account radiation in the first-order approximation amounts to solving the linearized regular 
Riemann-Hilbert problem (or the jump problem), 







m) 



m\e 



-2ikx 







r-\k), 



J 



where we have used the notation |6) = b = {bi, b^)'^. Solution of the above jump problem 
is obtained by integration and using the normalization condition 0±(A;) — * / as A; ^ oo. We 
obtain 



1 7 di 







V 



(6(£)|e-2^^" 

The contribution from radiation to the solution follows from the formula 
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oo 

(rad) _ _2 jj^ k(j)iN{k) = — f dk\ r{k) 



' e^^^\h{k))^ 



r-^(A:) 



As in derivation of the equations for the sohton parameters it is convenient to use the 
transformation with the matrix S. Using this transformation and formula ([75|) one can 
easily simplify the formula for radiation contribution. We get (ii'^'^^ = g;"'' + where the 
"parallel" and "perpendicular" parts of radiation are defined as follows 

qiW) ^ n_ / ^^|(A + ,tanhz)2e^^^(^|^(A)) + sech2ze-^-(^(A)|^^)}, (92a) 



gj^) = ^- / ^e^^\\ + ii^n\iz){giX)-e,{e\g{X))). (92b) 



— oo 



Here, for convenience, we have introduced the modified spectral parameter by setting 
k = rj\ + ^ and the modified spectral densities of radiation: 



9i 



(\) =e^'^'^^+^Hl{r|X + 0■ (93) 



The separation of radiation into parallel and perpendicular parts is due the following facts. 
While the perpendicular part ql^'^ satisfies the orthogonality property 

f:%(^) = 0, (94) 

1=1 

the parallel part of the radiation is given by the same formula as the radiation of the scalar 
soliton (multiplied by 9i), but for the averaged spectral density {9\g{X)). 

Consider the evolution equation (^) for the spectral densities of radiation. In the first- 
order approximation the term with T^n can be neglected. When deriving evolution equa- 
tions for the modified spectral densities in the first-order approximation one must take into 
account only the fast (or "integrable" ) evolution of the soliton parameters involved in the 
definition of g{X)- In particular, the modified spectral parameter A is t-independent. The 
only parameter which has the fast t-dependence in ( P^ is x. By differentiation of ( P5| ) and 
using (|7lD and (^) we obtain 
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where 



ririX) = 2[Xlm{u{h)} + Re{uj{h)} - uj{riX + 



(95) 



(96) 



The last two terms in equation ( pS]) come from Tin (0) (the second one is the contribution 
of (0)). They read 



T(II)(A) 



2r]{X^ + l) 



dz e 



-i\z 



(A — i tanh z) tq — sech z ro 



(97) 



Ti^)(A) = — / dze~^^"(A - itanh;z)(e-^(«/'')Vi - 6, 

2ri[X + i) J 



(98) 



Due to the definition of tq (^O]), the perpendicular component satisfies the identity 

E^^'^J^^ = 0- (99) 



1=1 



Integrating the equation for gi{X) with gi{X,t = 0) = and using the result in (p^ ) we 
arrive at the first-order correction to initially pure soliton solution. In this case we obtain: 



ITT 



J ^^^{(A + 2tanh^)V^^+^^^WS^"nA) 



where 



+ sech2;2e-^^^-^^'-W*7(ll)(A)}, 



g;^) = 6^^^ 7 -^(A + .tanh.)e'^^+^^^(^)S;^)(A), 

ITT J X- I 

— CO 



7^ 



1 



2r]{X^ + 1) 



dz e 



-iXz 



(A — ztanh2:)^fQ — sech^2;fQ 



(100a) 
(100b) 



1 r 

li ' = 2r^^x + i) J d^e-^''(A-^tanh;.)(e-^(«/'')^n-^4+)), 
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with 



(+) 



ro{z,T), '{z,t,X) = 



dTe-'^^^^^^ro{z,r), 











t 




dr e 








VI. CONCLUSIONS 



In construction of the perturbation theory our main idea is to use the Riemann-Hilbert 
problem associated with the integrable PDE for the nonhnear transformation of the perturbed 
PDE to the spectral space. The evolution equations for the spectral data follow from the 
evolution functional, an additional object one needs to introduce into the 1ST theory to 
account for perturbations. For a single vector soliton, the equations describing evolution 
of the soliton parameters and first-order radiation are given in explicit form. The method 
is not restricted to the first order only. For instance, the second-order equations can also 
be derived. The perturbation theory can be applied for description of dynamics of the 
spatial optical solitons, soliton pulses in the multispecies Bose-Einstein condensates, soliton 
propagation in optical fibre with the account of the arbitrary polarization of light pulses, 
and for many other applications of the multi-component soliton equations. 

In this paper we have restricted the consideration to the Zakharov-Shabat spectral prob- 
lem. However, the approach of this paper was successfully applied to other spectral problems 
as well [41-45]. There, the evolution functional was derived and the evolution equations for 
the spectral data were obtained. The overall result of this and the previous works on the 
perturbation theory based on the Riemann-Hilbert problem is that this approach always 
works. The explicit form of the evolution functional was the same for all considered spectral 
problems and, moreover, it undergoes only insignificant changes in the transition from the 
Cauchy problem to an initial-boundary value problem [45]. 
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APPENDIX A: SOME COMMENTS ON MULTIPLICITY OF ZEROS 

Here we explore in more detail the algebraic and geometric multiplicities of zeros of 
det $+(/;;) and det <I'I^(A;). For instance, consider the determinant 

det $+ = J+, A ... A J+„ A J_„+i A ... A J_^. (Al) 

None of the columns J±i{x, k) is equal to zero (otherwise, due to uniqueness of solution to 
equation (|l|) we would have J±i = for all x). Moreover, the minors in the wedge products 
of columns of J+ and J_ satisfy linear homogeneous equations, which follow from (|l|), section 



9^.J+i A ... A J+„ 



xJ-n+i A ... A J_^. 

Here M^^^ denotes a super-matrix, whose action on the wedge products of vector-columns 
is defined by the rule 

j 

M^^^^i A ... A = ^ ^1 A ... A M^p A ... A ^j. 

p=i 

Therefore, rankJ^^^ A ... A J+„ = rankci A ... A e„ = n, and rankJ_„_|_^ A ... A J_jv = 
ranke^+i A ... A eAr_„ = N — n. Hence, the only possibility for det ^+{k) = is the linear 
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dependence of the columns in J-^Hi and J-H2, e.g., at least one of the columns of J-H2 is 
given as a linear combination of the columns of J+Hi. 

Let Uj and dj be the algebraic and geometric multiplicities of the j-th zero kj, i.e., 

det $+(A;) = (A; - kjY^ipik), ip{kj) ^0, dj = N - Tank^+{kj). 

Writing down the Taylor expansion of ^+{k) about k = kj we immediately conclude that 

rank$+(A;j) > N - uj. (A2) 

Therefore, in general, the algebraic multiplicity is greater than the geometric one; trivially, 
they coincide for simple zeros. Representation (^) gives 

rank$_|_(A;j) > max(n, N — n). 

Hence the geometric multiplicity satisfies 

dj < N - max{n, N -n). (A3) 

We consider only zeros whose algebraic multiplicity is equal to the geometric one. In 
particular, if n = — 1, there can be only one vector in the null space of $+, i.e., dj = 1. 
Hence, in this case, zeros of det$+(A;) must be simple to satisfy the equal multiplicity 
condition. This condition requires that the only case of a multiple zero, det^+{kj) = 0, of 
order Uj is that there are precisely uj columns of J^{kj)H2 and k)Hi which are linear 
combinations of the columns of J^{kj)Hi and J^{kj)H2, respectively. 

The equal multiplicity condition can be guaranteed by the following constraint (imposed 
for some (x, t)) 

rank U+, \ = N. (A4) 

Indeed, for algebraic multiplicity uj, by the Taylor expansion, ( |A4| ) gives rank$+(/cj) = 
A^ — Uj. Conversely, if rank$+(A;j) = N — uj, for algebraic multiplicity uj, then there are at 
least Uj columns in d^+{kj)/dk independent from the columns of ^j^{kj). Hence ( |A4| ) also 
holds. Similar results are valid for multiplicities of zeros kj, j = 1, s, of det $I^(fc). 
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APPENDIX B: PROPERTIES OF THE REGULARIZATION MATRIX 



Here we derive the regularization matrix r(A;) and prove its properties (in the main, we 
follow Refs. 13 and 58). Dependence on the co-ordinates x and t is not important for this 
purpose and omitted. Consider one pair of zeros, say, kg and kg of det ^+{k) and det ^Z^{k), 
respectively. Let the vectors \p\^^) and {Pi^^^l, I = 1,...,//^, satisfying 

$4MbS^)) = o, (pp)|$:i(fc,) = o, (Bi) 

span the respective null spaces. Construct the following rational matrix functions 
where 

l,m=l 

and Ps is a projector: = Ps, rankP^ = z/^. It is easy to verify that XsW is inverse to 
Xs{k): Xs{k)Xs{k) = I- 

The determinant of Xs{k) is easily computed in some appropriate basis, where the pro- 
jector is represented by a diagonal matrix with Ug ones and N — Ug zeros on the diagonal. 
We get 



det Xs{k) = 




Hence, with such rational matrices we can factor out the s-th pair of zeros. Indeed, con- 
sider the products $+(A;)x7^(A;) and Xs{k)^Z^{k). These matrix functions are holomorphic 
in the upper and lower half planes, respectively (the poles are removable due to the identities 
( pID ). On the other hand, the determinants are non-zero for k = kg and k = kg, respectively; 
thus one pair of zeros is factored out. By introducing a sequence of such matrices, 

xAk) = I-^^Pv Xjik)=I + ^^P,, j = l,...,., (B2) 
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where Xj(^) = Xj ^{k), we factor out all zeros using the regularization matrix T and its 
inverse, where 

T{k)=xiik)x2ik)-...-Xsik). (B3) 
The projector Pj is given by the following formula 

Here the vectors lep"*) and (ep^| are related to the basis vectors of the null spaces of ^+{kj) 
and by triangular equations (if scanned starting from s down to 1): 

bP) = x;\k,) ■ x:\{k,) ■ ... ■ x7i(fc,)|eP), (B4) 

{e?^\XHi(kj) ■ xMkj) ■ ... ■ Xs(k,) = (pFY (B5) 

Due to Pj = Pj, these vectors satisfy the identities: 

xAkj)\e?) = 0, {ei^\xj\kj) = 0, l = l,...,u,. (B6) 

The regularization matrix T{k) can be made parameterized entirely by the vectors from 
the null spaces. Indeed, let us decompose T{k) and the inverse matrix into the partial 
fractions: 

j=l f^j j=l 1^ 1^] 

where due to ( p^ and ( pSj) we have 

B,=Y.\P?){-W ^,=f:|l;?)(#|. (B8) 

1=1 1=1 

From (|B7| ) and the identity FF"^ = F^^F = / it follows that 

F(A;,)bP) = 0, (pp)|F-i(I,) = 0, l = l,...,u^, J = l,...,s. (B9) 

These are the equations defining the unknown vectors \v\''^) and (fp''|. Indeed, rewriting 
59f) we have 
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i=l f^j m=l i=l '^j '^i m=l 

Inversion of these formulae gives 

i^i") =it ip2> . <"i"i =it -'),,,,„ m (BIO) 

j=l m=l i=l m=l 

Here the matrix D is defined by 




'i 



Substitution of (piCI| ) into ( pTf ) produces the needed formulae (p9|). 
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